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ABSTRACT 



In its motion response to a disturbance, a ship can be closely ap- 
proximated as a second order system. The solution of these second order 
equations is used to determine the motion stability of the ship. By ap- 
plying a step disturbance to a ship model and measuring the response, 
the linearity can be used to find the response to any disturbance. 

In this work a step disturbance was applied to a Mariner class 
hullform and the response measured. Prom this response the response to 
regular sinusoidal excitations was computed. The response to the sinu- 
soidal inputs was then used to determine the coefficients of virtual 
inertia and damping for the Mariner in pitch and heave. 

The derived results at zero speed were compared with values com- 
puted by strip theory, 
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CHAPTER I 



INTRODUCTION 

The equations of ship motion developed by Abkowitz (1) are a 
linearized set of second order partial differential equations. In the 
vertical plane the equations are 

(m-Z^.)z - - Z^z - (Z.^+mX^)e - (Z.+u^Z»)9 - (Z^+u^Z^)© = Z(ex) (l) 

(I - (M-+U M )e - (M^+u M. )e - (M.-+mK„)z - M-Z - M z = M(ex) (2) 

^ y 0 -' ^ e o z' 6 o z^^ '■ z G Z z ' ' ^ ^ 



z “ displacement in the vertical plane 
z = velocity in the vertical plane 
z “ acceleration in the vertical plane 
0 * angular displacement in pitch 

0 ® angular velocity in pitch 
'0* = angular acceleration in pitch 
Uq ® forward velocity of the ship 
m “ mass of the ship 

Z-- “ added mass coefficient 

z 

Z£ •= damping Coefficient 

Z^ ® static restoring force coefficient 

«= distance from origin to the center of gravity( equal to zero) 

G 

Z- " coupled virtual mass 
0 

Z* = rotary force derivative coefficient 
9 

Z “ static force derivative coefficient 

8 

I = moment of inertia 

y 

M*- “ added moment of inertia coefficient 

0 

M* = damping coefficient 
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Mq *> static restoring moment coefficient 

M-* •» coupled virtual moment of inertia coefficient 

M. ■ rotary moment derivative coefficient 

M “ static moment derivative coefficient 

z 

These equations show a coupling effect between pitch and heave. 
However, by decoupling these equations and neglecting second order 
effects, these equations can be reduced to 



The equations are now in a fom that is readily solvable. The 
solution of these uncoupled equations is highly significant in deter- 
mining the motion characteristics of a ship. These characteristics are 
Important in determining stability of motion as well as factors affect- 
ing the comfort and safety of personnel embarked on the ship. Before 
these equations can be solved, however, it is first necessary to deter- 
mine the hydrodynamic coefficients of added inertia and damping. 

The added inertia is a term that results from the fact that a ship 
moves in a dense medium and therefore must displace a heavy mass of 
liquid whenever it moves in any plane. This movement of the mass of 
liquid has the effect of essentially increasing the mass of the ship. 

The damping is a term associated with energy dissipation and not vis- 
cous effects. As a surface ship operates at the interface between water 



(m - Z-*)"z - Z.z - Z z = Z(ex) 
z z z 



( 3 ) 



(ly - M^.)e - M^e - Mg6 = M(ex) 
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and air, any vertical motion results in the creation of waves. The 
energy lost in wave creation acts as a damping force to reduce the os- 
cillations of the ship. 

Many methods, both theoretical and experimental, have been de- 
veloped to determine these coefficients. One such method was developed 
by Narita (11) at the Massachusetts Institute of Technology in 1962, 

It consisted of applying a known moment step disturbance to a hull form 
and measuring the transient response in pitch, while restraining the 
model from motion in heave. The resultant motion of the model was a 
transient wave of a decaying sinusoid form leading to the step. When 
this transient response was obtained, it was used to compute the theo- 
retical response of the model to any known input. This makes use of a 
property of a linear system, that when its response to a step or im- 
pulse input is known, its response to any other input may be found by 
use of the convolution integral. As a ship is closely approximated by a 
linear second order system, this principle is applicable. 

This property can be used to apply known regular sinusoidal inputs 
and determine the ship's response. This in essence produces the same 
results as could be obtained by making forced oscillation tests. The 
output of the regular sinusoidal inputs can then be used to calculate 
added mass and damping, 

Narita (11) computed added inertia and damping for the pitch 
motion only. In this work the procedure was extended to measure the co- 
efficients for pitch and heave. Two series of tests were run, one in 
which the model could pitch with heave restrained and one in which the 



modd. coiild heave with pitch restrained, 



The authors know of no previous tests conducted on the Mariner 
hull form for the purpose of determining coefficients of damping and 
added mass. However, a theoretical prediction for the zero-speed con- 
dition may be approximated. Porter (lU) has reported that the conformal 
transfomation of the unit circle represented by the equation 



where z is the plane of the actual ship section and C is the plane of 
the unit circle, can represent a, given station of the ship, N is se- 
lected and the coefficients a 2 nti are found which give the best fit to 
the actual section, Lewis (10) used this transfom with N “ 1 and 
Landweber and Macagno (9) extended it to N ® 2, In the present work, 
a method was used which can give a least mean square error approximation 
to the actual ship section if the number of coefficients is selected. 

The least number of coefficients which gave a reasonable fit for each 
section was used. 

After the coefficients of the transform were determined. Porter's 
method was used to obtain coefficients of damping and added mass for 
each station at several frequencies and the results were numerically in- 
tegrated throughout the length of the model by using strip theory to 
obtain added mass and damping for the raodd.. 
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CHAPTER II 



PROCEDUBE 

2,1 Model Characteristics 

The tests were conducted using a model of the Mariner class hull 
form, which has the characteristics shown in Table I, 



liable ^ 

CHARACTERISTICS OF THE MODEL 



Length between perpendicvilars 

Beam 

Draft 

Displacement in fresh water 
Block coefficient 
Radius of gyration 
Static restoring moment (pitch) 
Static restoring force (heave) 
Material of model 



5,28 feet 
,76 feet 
.27 feet 
Ii 2 ,l 4 . pounds 
.6125 

1.26 feet 

158,81| foot - pounds/radian 
lU8,80 pounds/foot 
Fiberglass, reinforced plastic 



2,2 Test Conditions 

The tests were conducted at the Massachusetts Institute of Tech- 
nology ship model towing tank. The tests were divided into two series, 
pitch with no heave and heave with no pitch. Each series consisted of 
tests conducted at Froude numbers of 0,0, 0,1, 0,2, 0,3. Further, the 
tests conducted at zero Froude number were conducted parallel to and 
perpendicular to the centerline of the towing tank in order to investigate 
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the effects of wall reflections. In the pitch - no heave, configuration, 
the ballast weights were shifted in order to determine the effect of 
changing the radius of gyration on the computed coefficients, 

2,3 Step Disturbance Application 

2.3.1 General 

As the tests were to be conducted in a liquid medium, it was de- 
cided that the high potential devices, as recommended by Narita (11), 
posed a danger to personnel working with the equipment and to other 
electronic equipment in the event of a wiring mistake or a faulty lead. 
For this reason low current and voltage devices were investigated, and 
it was decided that the best device might be a small electromagnet. The 
magnet chosen for these tests was a simple iron core wound magnet. For 
the purposes of these tests, the magnet required less than two amperes 
of current and less than two volts of potential. These values were not 
only safe to personnel in the event of a casualty, but were easily ob- 
tainable with standard direct current power supplies. Initially bat- 
teries were employed, but as the tests took a long time, because of de- 
lays in waiting for the tank to settle to a dead calm, the batteries 
rapidly lost their potential. In place of the batteries a Sorenson 
Model QR 18-1, ^A direct current power supply was used, 

2.3.2 Pitch - No Heave Configuration 

To apply the pitch disturbance, the model was mounted as shown in 
Figure I. A small metal plate was moxmted on the bow and on the stern 
of the mod^, A millimeter bar was used to mount the magnet. This 
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Figure I 

PITCH-NO HEAVE CONFIGURATION 
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enabled the operators to have greater ease in setting the disturbance at 
very close to one and a half degrees. As mentioned above, ’’the current 
was supplied to the magnet from a Sorenson Model QR 18-1, direct cur- 
rent power supply. The electric wiring arrangements were as shown in 
Figure II, By releasing the current to the magnet, the ship will let go 
under a step-like displacement, following a decasdng sinusoid coming to 
rest at a zero angle of pitch. 

As the exact moment is desired at which the ship actually reacts 
to the disturbance, a second circuit was constructed as shown in 
Figure m. When the ship begins to react, the magnet parts from the 
steel plate and produces a step reproduction on the recorder. The be- 
ginning point of this step can be used to detemine the starting point 
of the mode's response and thus increase the accuracy of calculation 
based on the response of the model, 

2,3»3 Heave — No Pitch Configuration 

In constructing the heave arrangements, it was necessary to devise 
a method of lifting the entire ship. It was also desirable to not in- 
troduce any forces tending to cause pitch, which would thus put undue 
strain on the heave rod. For this reason it was decided to mount a 
small metal plug directly on the heave rod and to lift the model by 
connecting the magnet to the metal plug. As the heave rod is moimted 
approximately at the longitudinal center of floatation of the model, es- 
sentially no pitch moment is introduced. The magnet was fastened to the 
carriage by two aluminum rods as shown in Figure IV, The circuitry is 
exactly the same as for the pitch configuration. 
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Figure II 

ELECTRIC WIRING ARRANGEMENT 
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Figure IV 

HEAVE-NO PITCH CONFIGURATION 
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lii both heave and pitch testing the magnet and metal connecting 
plate wejce mounted rigidly, so there was no concern for catching a fall- 
ing weight* Then to reconnect, it was only necessary to close the 
switch and lift the model so th^t the plate and the magnet came in con- 
tact, at which time the magnet would support the weight of the model. 
Then it was only necessary to wait until the tank settled to an absolute 
calm before beginning the next run, 

2,li Measurement of the Step Response 
2,U,1 Pitch 

In order to increase the validity of linearity of motion theory, 
an initial displacement of only one and a half degrees was used. To get 
the maximum accuracy in measuring this small variation, two methods were 
investigated, use of a linears 3 m and use of a pitch bearing. To use a 
linearsyn, maximum excursion of motion is desired. This would be at the 
ends of the model. However, the motion here is circular enough to cause 
non-linear motion of the linearsyn and create inaccuracies. Also, the 
vibrations near the ends of the model were found to be the greatest, and 
this too would have an adverse effect on measurements made by a linear- 
syn, 

A pitch bearing located right at the longitudinal center of 
floatation has the advantage that it can measure angular motion right 
at the point about which pitch takes place. The pitch beai*ing is a 
rotary motion displacement pickoff of a differential transfomner type 
and is highly sensitive to any rotary displacement. Therefore, it was 
decided to employ the pitch bearing for the pitch measurements. 
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2,h»2 Heave 



To measure heave, the linearsyn was perfectly suited. It is a 
linear motion displacement pickoff of the differential transformer type. 
It consists of a coil assembly and movable magnetic core and is capable 
of detecting motions as small as 0.000001 of an inch between the core and 
coil assembly. Mounted directly onto the heave rod assembly, the 
linearsyn gives an accurate reading as to the movement of the model in 

■s 

the heave plane. 

To increase the validity of model linearity in heave motion, a 
disturbance of only foujvtenths of an inch was applied to the model. 

This proved sufficient to obtain several well-defined cycles of oscil- 
lation® 

2,5 Test Procedure 

After the model was mounted on the carriage and circuitry arranged, 
the magnet was connected to either displace the model in pitch or heave, 
depending on which tests were being conducted. Tests were then carried 
out in two groups, pitch and heave, at Froude Numbers of 0,0, 0,1, 0,2, 

0,3. 

With the model displaced it was necessary to wait until the tank 

I 

had settled to an absolute calm before beginning a test. Then the model 
was brought up to speed and the switch opened, stopping power to the 
magnet. At that instant the ship would fall toward zero in a decaying 
sinusoidal motion, characteristic of a second order system. This motion, 
as well as the exact instant at which the magnet parted from the metal 
plate, was recorded on both the Sanborn recording paper and on magnetic 
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tape for later analysis. To begin the next run, the model had only to 
be towed back to the opposite end of the tank, the magnet re-engaged 
with the metal plate, and then allowed to remain suspended until the 
tank settled. 

For the zero speed runs the model was placed at the center of the 
tank, first parallel to and then perpendicular to the centerline. Here 
the procedure was the same as before except that there was no forward 
motion of the model. 
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CHAPTER III 



ANALYSIS 

3.1 Background 

The analysis procedure used In this work was essentially that used 
by Keivin and Narita (5)* A step input is applied to a system, and re- 
sponse is represented by 

P(t) = i (1 + R(t)) (6) 

where C is either the static restoring force in heave or the static re- 
storing moment in pitch, R(t) represents the transient portion of the 
response and l/C is the final value the response seeks after the tran- 
sient has died out. 

To get the response to a regular sinusoidal input, the convolution 
integral is used. 

00 

Vo(t) = f V.(t) H(t - t) dT (T) 

H(t) represents the impulse response of the system. However, as the step 

response, S(t), is knovm, it may be used with the input in the form 

dV. (t) 

-if- 

and the same resvilts will be obtained. From this, the response of the 
model to sinusoidal inputs in either heave or pitch may be obtained. 

To obtain the damping and -irirtual mass coefficients, we represent 



our second order equations as foUowss 



H ^ ®P if * V' “o 

at 

^ ^ ^ ^ S' 



(9) 

( 10 ) 



The solution of these equations will be of the fom 

a sinout + b cosyjt 



from which 






Cp - M^a^/(a/ . t/) 



0) 



^p = 



-Vp 



w(ap^ + bp^) 



^ = 



- VlA^ " 



03 



= 



-Z b„ 
o H 



/ 2 ^2. 

®(a„ + bj, ) 



( 11 ) 

( 12 ) 

(13) 

(1^) 

(15) 



3.2 Static Coefficients 

Equation (6) indicates that the step response tends to l/c as the 
transient dies out. In the initial simplification of the equations^ it 
was assumed that this was a static term and not speed dependent, ■ Ger» 
ritsma (2) shows that this term is speed dependent and that this depends 
ency is further a function of the hull form. While the variation in 



speed is not great at very low Froude numbers 5 it does approach approxi- 
mately 10 percent at a Froude number of three-tenths , This is still not a 
large effect, but it may exert a measurable influence on the values of 
added mass and damping calculated at the various Froude ntimbers. 

Accordingly, tests were conducted to determine the change in the 
static restoring coefficients with speed. In pitch a known moment was 
applied to the ship and the pitch angle was measured. As the pitch 
velocity and acceleration were 'both zero, the static restoring coeffi- 
cient could be found by dividing the moment by the resultant angle. The 
static restoring moment at zero speed was computed to be 158, 8U ft, -lb,/ 
radian. This same procedure was repeated at Froude numbers of 0,1, 0,2, 
and 0,3, and a variation of less than 1 percent was observed at each 
speed. Therefore, the static coefficient at zero speed was used for the 
analysis at each of the Froude numbers, 

A known heaving force was then applied to the model, and heave was 
measured at Froude numbers of 0,0, 0,1, 0,2, 0,3. In this case the re- 
sults followed the form of those found by Gerritsma as shown in Table II, 



Table H 



STATIC RESTORING COEFFICIENTS IN HEAVE 



Froude Number 



2 ^ 



,1 



0 



2 



Hi8,8 Ibs/ft 
lla,0 Ibs/ft 
lii2,8 Ibs/ft 
139.2 Ibs/ft 
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The variation is not large, but is significant enough to warrant 
using the value found at its corresponding Froude number* 

3*3 Digitizing Rate 

The data as obtained from the experiments was recorded in analog 
form on magnetic tape* As the convolution integral was performed by 
series methods, it was necessary to have the data in digital form* Pre- 
vious experiments conducted by Kerwin and Narita (5) showed the results 

■5 

of added mass and damping to be highly dependent onthe digitizing rate 
used, with damping being the more critical* From their tests of a 
known second order system, it was found that a digitizing rate of one 
kilocycle per second gave exact results for added mass and to within 
three-tenths of 1 per cent for damping. 

Accordingly, a digitizing rate of one-thousandth of a second was 
used. To digitize the data at this rate, it was necessary to record a 
one kilocycle per second square wave on the tape, as the process used 
with the IBM 709h computer reads the value of data everytime there is a 
positive to negative vertical change in the digitizing signal* 

As mentioned earlier, the response of the model was recorded as 
well as a signal that showed the instant the magnet parted. These sig- 
nals were recorded first, and the square wave was recorded afterward. 

The square wave was recorded insuring that it started prior to the in- 
stant of the magnet parting signal. By this procedure and obtaining the 
digitized data for both the response data channel and the parting signal, 
the parting signal data could be scanned, and the exact instant of mag- 
net parting could be found. On the other channel, this was the starting 
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point for which the response signal data was used in the analysis pro- 



gram. 



As mentioned above, to get maximum accuracy, a digitizing rate of 
a thousand points per second was used to obtain the best representation 
of the model response data, Tn convolving this data with the input 
sinusoid, a series method was used. This involved shifting the sinu- 
soid in reference to the response data. The response to the sinusoid 
was computed as a series, and tjie final value summed. In shifting the 
waveform, it was found that the fine spacing of .001 seconds was not 
required. The analysis program was altered to determine what interval 
could be used and still produce the same results that one-thousandth of 
a second piDduced, The main purpose of this alteration was to speed 
up the analysis program and reduce computer time used in obtaining re- 
sults, Table III shows the results of increasing the convolution 
interval at several frequencies or corresponding periods of excitation. 
These values were obtained using ideal data for the second order 



From Table HI it can be seen that increasing the convolution in- 
terval actually improved the accuracy of results at all except the high- 
est frequency, period equal to 0,6Ii seconds. In this case the added in- 
ertia stayed approximately the same, but the damping became worse by 
seven- thousandths. The over-all results did show that larger convolution 
intervals gave accurate results. Accordingly, larger intervals were used 



system 




(16) 



dt 
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in the analysis, depending on the input excitation frequency 

Table in 

EFFECT OF CONVOLUTION INTERVAL ON COEFFICIENTS 




Period 


Convolution Interval 


Added Inertia 


Damping _ 


(Seconds) 


(Seconds) 


(Slug Feet^) 


(Slue Feet^) 
Second 


1,28 


,001 


2.50.7 


7.581 


1,28 


,010 


2.503 


7.550 


1,28 


.010 


2.503 


7 . 5 U 9 


1,28 


,020 


2.503 


7 . 5 U 9 


1,28 


.032' 


2.502 


7 . 5 U 6 


1,28 


0 

0 

• 


2.502 


7 . 5 U 5 


2,56 


,001 


2.511 


7.658 


2,56 


,010 


2.506 


7.602 


2.56 


.016 


2.506 


7.601 


2,56 


.020 


2.506 


7.600 


. 6 U 


,001 


2 .ii 98 


7.529 


.6it 


.010 


2.U99 


7.536 


.6U 


,020 


2 . 1^99 


7.536 



3 «li Theoretical Analysis 

In the present problem offsets for 21 stations of the Mariner hull 
were read off a ■5-inch-scale body plan, and by iteration on the IBM 709U 
computer, the coefficients »2n+l which gave the best fit to the actual 
sections were determined. No more than five such coefficients for any 
one station were used because the ph5rsical significance of coefficients 
beyond becomes unclear. The values of the coefficients giving the 
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best fit for each section are tabulated in Table IV, 

After the coefficients «2n+l determined, they were used as 

inputs to a computer program which utilizes Porter’s method as explained 
in reference ih and calculates coefficients of damping and added mass 
for each station. The results for each station were integrated over the 
length of the model by using strip theory to give the coefficients for 
the model. 
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1 

2 

3 

h 

5 

6 

7 

8 

9 

10 

11 

12 

13 

Ik 

15 

16 

17 

18 

19 

20 



Table TV 



SECTION TRANSFORMATION COEFFICIENTS 



Half-Beam 


Draft 


Section Area 
Coefficient 


0,20 


27.0 


7. ii 07 


h ,30 ■ 


27.0 


0.999 


9.75 


27.0 


0.818 


16.15 


27.0 


0.779 


22,70 


27.0 


0.791 


28.50' 


27.0 


0.830 


33.25 


27.0 


0.86 U 


36.30 


27.0 


0.9 lii 


37.70 


27.0 


0.958 


38.00 


27.0 


0.982 


38,00 


27.0 


0.988 


38.00 


27.0 


0.979 


38.00 


27.0 


0.965 


38.00 


27.0 


0.932 


37.75 


27.0 


0.868 


36. ii 5 


27.0 


0,787 


33 .hO 


27.0 


0.700 


28.25 


27.0 


0,609 


21.05 


27.0 


0.519 


12,20 


27.0 


0,390 


1.80 


0.85 


0.020 
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Table IV ( continued) 



Station 




!3 






0 


-. 8^755187 


-.11659060 


-.Ol 89 l 6 ii 9 U +.0065186719 


-.00037253161* 


1 


-. 66039 U 


-.102210 


-,006052 +,02Hli2 




2 


-.h 739 ia 


-.Oii 8562 . 


♦.012796 +.0172^6 


+. 0061*57 


3 


-.273985 


+.O2U895 


+. 02^931 +.015372 




h 


-. 112^73 


-.017126 


♦.019732 +. 007 U 28 


+.007161 


5 


+.001188 


-.032628 


♦.019083 +.000588 


+.005890 


6 


+.076083 


-.053I492 


+.016873 +.000628 


♦ .005291* 


7 


+.112537 


-.078719 


+.010882 -.001901 


+.001656 


8 


+. 1^3558 


-.103009 


♦ , 00 i ;7 8 h 




9 


+.150UIO 


-.118791 


-.001681 -.002359 




10 


+,150708 


-.125100 


-.003798 -.002210 


+.000777 


11 


+.151138 


-.123101 


-.OO32U7 -.003000 




12 


+.l 508 h 3 


-.110398 


-.001079 -. 001^638 




13 


+.li; 85 oo 


-.085690 


+ .005^38 -.OOli 678 




lU 


+ .lli 0856 


-. 0 h 8 ii 39 


+.015873 -.OOla 78 


+.000558 


15 


+. 12 U 07 h 


-.002855 


+.023ii6l -.001276 


+.000786 


16 


+. 08 lli 37 


+ ,Oii 96 lO 


+.029596 -.001735 




17 


-.022515 


+,101750 


+. 01*2262 -. 0001*02 


+.005171 


18 


-.220882 


+.138063 


+.078202 +.001977 


+.001509 


19 


-. 551 i 295 


+.110698 


+ .101173 +.035001* 


+.020560 


20 


♦.319277 


-, 08 Ii 2 l 5 


+ .051*789 -.052018 


♦.019595 
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CHAPTER IV 



RESULTS 

i;.l Heave-No Pitch Condition 

To obtain values of added mass and damping in heave, four tests 
were made at each of the Froude numbers desired. A typical response 
curve is shown in Appendix B , These curves are presented to show 
wall reflection effects, but they also show the typical response at zero 
speed. At forward speeds greater than zero the model response curves 
were much the same, but considerable vibration was also introduced be- 
cause of the towing carriage oscillations, 

1;.1,1 Virtual Mass 

The tests conducted at each speed were averaged to obtain the vir- 
tual mass. This average vklue was then nondimensionalized by dividing 
it by the mass of the model. Each of these values is shown in Table V 
as a function of excitation frequency, Furthemore, each of the ob- 
served values was compared with the average to obtain the worst excur- 
sion, This worst excursion from the average is expressed as a percent 
of the average in Table V, 

In most cases all points were veiy close together, and percentage 
errors were small. However, occasionally a derived value was at vari- 
ance with the rest of the results and its error appeared quite large. 
This took place at those values farthest from the model’s natural fre- 
quency of approximately 7,6 radians per second. 
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The symbols used in Table V are as follows; 

w radian frequency of excitation 
A virtual mass 
A’ nondimensional added mass 
g gravity 

A displacement of the model 

The average values of the virtual mass are plotted in Figure V 
for each of the four Froude numbers. Figure VI is a comparison of the 
zero speed values obtained with the model parallel to and perpendicular 
to the centerline of the tank, 

li,l,2 Damping 

The same procedure was used with damping as was used with added 
mass. The results of the average values^ as well as the maximum percent 
error, are shoim in Table VI. The damping values were nondimens ional- 
ized by multiplying the damping obtained by the square root of the 
length of the model times gravity divided by the displacement of the 
model. 

Errors in damping were very small, but were generally slightly 
higher than for added mass. However, near the natural frequency of the 
model, the errors were very small. Increasing as the frequency was in- 
creased or decreased. 

Additional symbols introduced in Table VI are; 

B damping 

B’ nondimensionalized damping 

L length of model 



The average values of damping are plotted in Figure VII for each of 
the four Froude numbers. Figure Vm is a comparison of the zero speed 
values obtained with the model parallel to and perpendicular to the cen- 
terline of the tank, 

U.2 Pitch-Mo Heave Condition 

To obtain the values of the added inertia and damping coefficients 
in pitch, six tests were conducted at each of the Froude numbers. Three 

ft 

tests were conducted with the bow initially up, and three tests were con- 
ducted with the stem initially up. The response curves were the same as 
for the case of heave, in that each was a decaying sinusoid. At speeds 
other than zero there was a large amount of vibration in the response, 
paralleling the situation found in heave, 

li,2«l Virtual Inertia 

The same analysis program was used to compute the value of virtual 
inertia in pitch as was used to compute the virtual mass in heave. The 
only changes were to change several constants and the nondimensionalizing 
factor. Virtual inertia was nondimensionalized by multiplying the value 
found by the value of gravity divided by the displacement times the 
length squared. The values of virtual inertia found at each of the four 
Froude numbers are shown in Table VH, In this table A designates the 
value of virtual inertia. Figure IX is a plot of the virtual inertia 
versus excitation frequency for each of the four Froude numbers. Fig- 
ure X is a comparison of the zero speed condition with the model paral- 
Id. to and perpendicular to the centerline of the tank. This figure 
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also shows the effect on the virtual inertia of changing the radius of 
gyration of the model from 1,26 feet to 1,17 feet, 

U.2,2 Damping 

The damping was computed using the same analysis program as for 
heave with the appropriate constants changed. It was rKondimensionalized 
by multiplying by the square root of the length times the value of grav- 
ity and dividing this square root by the displacement times the length 
squared. The values of damping as a function of Froude number and fre- 
quency are shown in Table V3^, The damping is plotted in Figure XI to 
show its variation with Froude nvttaber. Figure XU is a plot of the zero 
speed conditions, also showing the effect of the reduced radius of gyra- 
tion, 

i;,3 Theoretical Results 

Based on the coefficients of transfomation as presented in Table IV, 
Porter's theory was applied to determine the coefficients of damping and 
virtual mass in both pitch and heave. These coefficients were determined 
at nine frequencies for each station, and the results were numerically 
integrated by strip theory to obtain the coefficients of virtual mass and 
damping for the modd.. These results are presented in Table IX, These 
same results are plotted in Figures XIIl-XVI to compare their values with 
the experimentally obtained zero speed values. 



- 26 - 











r 











Table V 

VIRTUAL MASS COEFFICIENT IN HEAVE 



(1) Froude Number 


0.3 






w 


A 


A* 


Max, ^ Error 


h,13h 


1.951 


1,1+81 


10,08 


U.909 


1.71+1 


1.322 


8.73 


6.01+2 


1.821+ 


1.386 


1.01+ 


7.11+0 


1.788 


1.358 


2.13 


7.851+ 


1.752 


1.330 


1.1+8 


8.727 


1.718 


1.303 


.93 


9.817 


1.702 


1.292 


.65 


11.220 


1.759 


1.336 


1.1+2 


12.083 


1.81+0 


1.398 


2.1+1+ 


13.090 


1.891 


1.U37 


3,75 


(2) Froude Number 


0.2 






w 


A 


A* 


Max, $ Error 


1+.131+ 


3.01+5 


2.313 


ll+,86 


1+.909 


2,308 


1.751 


9.50 


6,01+2 


2.129 


1.617 


8,65 


7.11+0 


1.959 


1,1+86 


6.1+8 


7.85U 


1.838 


1.395 


5.33 


8.727 


1.732 


1.316 


1+.27 


9.817 


1.651 


1.253 


U.18 


11.220 


1.601 


1.217 


3.26 


12,083 


1,625 


1.231+ 


7.08 


13.090 


1.678 


1.273 


8,11 
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Table V ( continued) 



(3) Froude Number 


0.1 






w 


A 


A* 


Max, % Error 


k.l3h 


3 .U 76 


2.619 


1+.21 


h.909 


2.692 


2,0ia 


2,71 


6.0U2 


2.282 


1.731+ 


7.50 


7.1ii0 


2.113 


1.605 


5.35 


7.85U 


1.978 


1.501 


2.68 


8.727 


1.81+5 . 


1.1+01 


3.63 


9.817 


1.757 


1.332 


l+.io 


11.220 


1.711+ 


1.301 


5.51+ 


12.083 


1.751+ 


1.331 


5.70 


13.090 


1,862 


I.I+II+ 


3.1+1+ 


(li) Froude Number 


0,0 (Parallel) 






w 


A 


A‘ 


Max, % Error 


U.13ii 


3.876 


2.938 


.10 


lt.909 


2.81+5 


2.159 


1.23 


6,01+2 


2.273 


1.723 


2.86 


7.lli0 


2,186 


1.659 


2,06 


7.851+ 


2.051+ 


1.555 


1,80 


8.727 


1.923 


1.1+59 


1.35 


9.817 


1.870 


1.1+18 


1.07 


11.220 


1.905 


1.1+1+1+ 


.68 


12,083 


1.879 


1.1+22 


.30 


13.090 


1.993 


1.511 


.71 
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Table V (continued) 



{$) Froude Number 0,0 ( Perpendicular) 



w 


A 


A» 


Max, % Error 


U.13U 


2.923 


2.220 


6.88 


h,909 


2.696 . 


2,01x2 


7.03 


6,Oh2 


2.267 


1.719 


ii.lO 


7.IU0 


2.097 


1,590 


3.39 


7 . 851 i 


2,030 


1.539 


li .28 


8,727 


1.958 ' 


l.Wli 


ii .75 


9.817 


1.900 


l.iiiiO 


5.21 


11.220 


1.901 


i.Ua 


6 . 81 , 


12,083 


1 . 9 li 9 


1 .U 78 


6,57 


13.090 


2.077 


1.572 


3.90 
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(1) Froude Nvroiber 
w 

U«13i; 

U.909 
6,0ii2 
7.1U0 
7.85U 
8.727 
9.817 
11.220 

12.083 

13.090 

(2) Froude 
w 

U.13U 
U.909 
6.0U2 
7.1U0 
7.85U 
8.727 
9.817 
11.220 

12.083 
13.090 

- 30 . 



Max . % Error 
3.67 
1.38 
U.62 
2.37 
3.92 
3.61 
3.12 
U.OO 
H.12 
' 9.27 



B* 


Max. % ] 


1.188 


10.78 


1.621 


15.U8 


1.509 


15.32 


1.828 


11,88 


1.910 


7,62 


1.931 


It.OO 


1.901 


8.27 


1.731 


U6.20 


1.578 


U9.70 


l.ii99 


100.00 



Table 71 

DAMPING COEFFICIENT IN HEAVE 



.3 



B 


B‘ 


7.333 


2.253 


6,180 


1.900 


5.991 


i.8ia 


6,l49it 


1 . 99 U 


6.637 ' 


2.039 


6.655 


2,Oit5 


6,ii20 


1.973 


5,929 


1,821 


6.173 


1.897 


7.295 


2.2ia 



Number .2 
B 

it,8U2 

5.279 

ii.910 

5.913 

6.a3 

6.281 

6.182 

5*636 

5.138 

U.87U 
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Table VI ( continued) 



(3) Froude Number ,1 



w 


B 


B« 


Max. % Error 


U.13U 


5.U15 


1.66ii 


17.67 


1.909 


7.302 


2.2U5 


lU.oo 


6.0U2 


6.U29 


1.976 


7.11 


7.H1O 


6.832 


2.100 


5.1U 


7.851i 


7.063 


2.172 


3.88 


8.727 


7.0U5 , 


2.163 


l.lU 


9.817 


6.i*63 


1.987 


1.89 


11.220 


5.2U8 


1.612 


9.91 


12.083 


ii.305 


1.322 


16,00 


13.090 


U.I493 


1.380 


27.21 


(ii) Froude Number 


0 (Parallel) 






w 


B 


B‘ 


Max. % Error 


U.13U 


5.911 


1.812 


6.96 


li.909 


8.1a9 


2.582 


5.51 


6.OU2 


6.660 


2.OU3 


2.12 


7.IU0 


7.126 


2.186 




7.85U 


7.390 


2.267 


.03 


8.727 


7w03U 


2.159 


.lU 


9.817 


6.112 


1.876 


.U9 


11.220 


5.550 


1.702 


.80 


12.083 


ii.991 


1.531 


.20 


13.090 


3.U3I 


1.052 


ii.20 
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Table VI ( continued) 
(5) Froude Number 0 ( Perpendicular) 



w 


B 




Max. % Error 


ii.l3li 


5.325 


1.635 


10.13 


h ,909 


7.1t30 


2.281 


10.52 


6.0U2 


8.1i6l 


2.600 


7.97 


7 .II 1 O 


8.30U 


2.551 


1.13 


7.85U 


8.2li2 


2.530 


1.72 


8.727 


8.005 


2.1i58 


1;.72 


9.817 


7.279 


2.236 


10,86 


11.220 


5.806 


1.782 


22,10 


12.083 


U.713 


IMS 


26.52 


13.090 


3.773 


1.158 


li5.30 



- 32 - 



I 




Table VII 



VIRTUAL INERTIA COEFFICIENT ni PITCH 
(1) Froude Number «3 



w 


A 


A’ 


Max, % Error 


U.13U 


2.I1IO 


,066 


7.3I4 


h,909 


2.501 


,068 


6,liO 


6,0U2 


2,286 


CM 

0 

• 


3.15 


7. Ho 


2,538 


,069 


1.03 


7.85h 


2.52I4 ' 


.069 


1.38 


8.727 


2.503 


,068 


1.68 


9.817 


2,622 


.071 


3.38 


11.220 


2.399 


.065 


3.71 


12.083 


2,593 


H 

0 

• 


U.32 


13.090 


2.1i20 


,066 


6.32 


(2) Froude Number 








w 


A 


A' 


Max. % Error 


U.13U 


2.986 


.081 


2.11 


I1.909 


2.528 


.069 


l.il2 


6.0lt2 


2,598 


.071 


2.38 


7. Ho 


2.667 


.073 


1.35 


7.8H 


2.605 


,071 


1.11 


8.727 


2.558 


.070 


.78 


9.817 


2.593 


.071 


,10 


11.220 


2 .li 3 U 


.066 


.70 


12.083 


2.55U 


.069 


I4.67 


13.090 


2.U03 


.065 


5.86 
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Table VII (continued) 



(3) Froude Number ,1 



w 


A 


A« 


Max, % Error 


ii.l3U 


U.020 


.109 


a.65 


U.909 


3.538 


.096 


3.7U 


6.0i;2 


3.019 


.082 


U.ii5 


7.11i0 


2.715 

-a 


.07li 


3.17 


7.85U 


2.551 


.069 


2.70 


8.727 


2.U12 ^ 


.066 


h,02 


9.817 


2.316 


.063 


6,09 


11.220 


2.200 


.060 


7.31 


12.083 


2.150 


.058 


9.35 


13.090 


2.181 


.059 


9.81 


( U) Froude Number 


0 (Parallel) 






w 


A 


A« 


Max, % Error 


ii.l3li 


3.623 


.099 


.ia 


li.909 


3.591 


CO 

ON 

0 


,h2 


6.0U2 


3.167 


.086 


.36 


7.11;0 


2.821 


.077 


.32 


7.85ii 


2.656 


.672 


.31 


8.727 


2.519 


.068 


.2U 


9.817 


2.ia8 


,066 


.17 


11.220 


2.279 


.062 


.32 


12.083 


2.199 


.060 


.10 


13.090 


2,li49 


.058 


.19 
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Table VII ( continued) 



(5) Froude Number 


0 (Perpendicxilar) 






w 


A 


A» 


Max. % Error 


U.13U 


1.999 


.05I4 


3.10 


U.909 


3.296 


.090 


i.Uo 


6.0U2 


3.017 


.082 


1.28 


7.1U0 


2,622 


.071 


1.37 


7.85U 


2.U80 


.067 


1.59 


8,727 


2.UO8 ^ 


.065 


1.7U 


9.817 


2.U9I 


.068 


2.69 


11.220 


2.273 


.062 


U.53 


12.083 


2.3li5 


.06U 


3.67 


13.090 


2.553 


.069 


1.61 


(6) Froude Number 


0 (Reduced R) 






w 


A 


A' 


Max, % Error 


i^.l3U 


3.079 


.08U 


.29 


h ,909 


3.062 


.083 


.^6 


6.OU2 


2.67U 


.073 


.10 


7.IU0 


2.357 


.06U 


,U7 


7.85ii 


2.235 


.061 


.76 


8.727 


2.I3U 


.058 


1.03 


9.817 


2.O7U 


.056 


1.30 


11.220 


2.037 


.055 


1.67 


12.083 


1.998 


.05U 


2.20 


13.090 


2.0li9 


.056 


2.30 
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Table VIII 



DAMPING COEFFICIENT IN PITCH 



(1) Froude Number ,3 



w 


B 


B' 


Max, % Error 




6.651 


.073 


32.30 


I1.909 


8.5lii 


.09U 


11.30 


6.0ii2 


6.157 


,068 


13.95 


T.liiO 


5.887 


.065 


8.92 


7.851i 


5.876 


.065 


ii.95 


8.727 


5.155 


.057 


2,26 


9.817 


5.130 


.057 


2.50 


11.220 


ii.3ii3 


00 

0 

• 


52.60 


12.083 


2.021 


,022 


100.00 


13.090 


ii.762 


.052 


iiO.50 


( 2 ) Fro ude Numbe r 








w 


B 


B' 


Max. % Error 


h , 13 h 


8,52ii 


,09ii 


U.ii2 


li.909 


7.82ii 


.086 


2.58 


6.0I42 


5.U68 


.060 


6,70 


7.11i0 


5.863 


.065 


3.11i 


7.85ii 


5.823 


.06U 


.86 


8.727 


5.189 


.057 


1.00 


9.817 


U.98I 


.055 


3,1*8 


11.220 


3.588 


.0ii3 


13.60 


12.083 


2,921 


.032 


21.60 


13.090 


2.003 


.022 


100.00 
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Table VHT ( continued) 



(3) Froude Number 


0.1 


w 


B 




3.087 


li.909 


1.316 


6 , 0 hZ 


5.35ii 


7.H+0 


6.0I4.8 


7.85I1 


6.318 


8.727 


^.900 . 


9.817 


S,hSh 


11.220 


ii.665 


12.083 


3.328 


13.090 


1.276 


(U) Fronde Number 


0 (Parallel) 


w 


B 


U.13ii 


.820 


ii.909 


3.99U 


6.0li2 


5.508 


7.11i0 


6.173 


7.85ii 


6.091 


8.727 


5.637 


9.817 


5.102 


11.220 


U.205 


12.083 


3.155 


13.090 


1 , 1 x 91 





Max. % Error 


.03li 


IU.95 


00 

• 


1,01 


.059 


1.36 


.067 


2.U8 


.070 


l, 2 h 


.06^ 


3.72 


.060 


5.81 


.051 


9.60 


.037 


13.36 


.oiU 


20.82 



B* 


Max, % Error 


.009 


.51 


.Oiili 


l.Oii 


.061 


,hZ 


.068 


.53 


.067 


.38 


.062 


.55 


.056 


2.66 


.0U6 


2.72 


.035 


2.ii0 


.016 


lh,hO 
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Table VIII (continued) 



(5) Froude Number 0 (Perpendlcvilar) 



w 


B 


B’ 

% 


Max. ^ Error 


Ii.l3l 


-2.723 


O 

o 

• 

1 


l.liO 


h.909 


1.351 


.015 


IU .36 


6.0U2 


6.197 


.068 


2.57 


7.11i0 


7.333 


.081 


l.Oii 


7.851t 


7.080 


.078 


.10 


8.727 


6,353 


.070 


.88 


9.817 


5.552 


.060 


1.15 


11.220 


7.1U5 


.079 


9.6U 


12.083 


3.53U 


.039 


10.83 


13.090 


10.385 


.111 


23.60 


( 6) Froude 


Nvmber 0 (Reduced R) 






w 


B 


B‘ 


Max. % Error 


b.l3U 


1.81i7 


.020 


16.83 


h.909 


3.886 


.013 


5.90 


6.012 


6.085 


.067 


2.26 


7.11t0 


6.550 


.072 


l.ii3 


7.85U 


6.1i55 


.071 


1.25 


8.727 


6.157 


.068 


.99 


9.817 


5.566 


.061 


.Ul 


11.220 


5.395 


.059 


1.26 


12.083 


1.701 


.052 


5.10 


13.090 


ii.578 


.050 


1.75 
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Table IX 



THEOHETICAL VIRTOAl MASS AND DAMPING 



Frequency 


Virtual 

Mass 


Damping 


Virtual 

Inertia 


Damping 


li*n 


2.823 


8.U31 


li.Oli; 


9.605 


^.82 


2.it2U 


8.5UO 


3.U68 


10.231 


7.13 


2.295 


7.703 


3.230 


9.857 


8.23 


2,257 


, 6.733 


3.139 


9.272 


9.20 


2.260 


5.933 


3.105 


8.536 


10.08 


2.281 


5.050 


3.073 


7.826 


10.89 


2.310 


U.153 


3.115 


6.378 


11.6U 


2,3liO 


3.893 


3.13U 


6.588 


12,80 


2.369 


3.086 


3.157 


5.U36 
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VIRTUAL INERTIA IN PITCH AS A FUNCTION OF FROUDE NUMBER 
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FIGURE XVI 



COMPARISON OF THEORETICAL AND EXPERIMENTAL 
ZERO SPEED DAMPINGT n PITCH 
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Chapter V 

DISCUSSION OF RESULTS 



5.1 Experimental Errors 

The resultant curves of virtual mass and damping as shown in Fig- 
ures V - XII show the respective values for the Mariner, In comparing 
these curves with those of other similar type large hullforms, these 
curves show the same trends in their values. These curves are gener- 
ally believed to give the values for the Mariner at the various speeds 
for the given excitations. However, there are several experimental 
errors or limitations which may limit their accuracy. 

Initially the magnet was supported by a wire. When the tests 
were conducted, the random vibrations were so great as to almost ob- 
scure the actual record of the model’s motion. For this reason stiff 
supporting members were used. This greatly reduced the vibrations, but 
the inherent vibrations in the towing carriage structure and members 
could not be eliminated. These vibrations became greater at the higher 
speeds. They were also affected by the rate at which the model was 
brought up to speed, being much greater when the model was quickly ac- 
celerated to the desired speed. In general the vibrations were high 
frequency with little effect on the results, but since the vibrations 
were so random, it is likely that there were enough low frequency com- 
ponents to affect the res\ilts. 

Another source of error was introduced by the procedure used in 
recording the results. Initially the response was recorded at 3.75 
feet per second. At this low speed a one-kilocycle-per-second square 
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wave would not reproduce on the tape* In transferring the response to 
a tape at 30 feet p>er second, it was observed that noise, primarily 
sixty cycles per second, was picked up. The sixty cycle noise should 
have had no effect on the results, but any lower frequency components 
would have. This effect was considered to be very minor. 

As mentioned earlier, a second circuit was constructed to record 
the actual instant at which the model began to drop, as evidenced by 
the magnet parting from the metal plate. By digitizing this response, 
it was hoped to be able to pinpoint the proper point on the model re- 
sponse data and not introduce any errors by starting the analysis at 
the wrong instant. The system was generally satisfactory, but it was 
sometimes difficult to be sure if the chosen starting point was within 
two or three thousandths of a second of the actual one. To find the 
effect of missing the proper starting point, some of the data was ana- 
lyzed at points away from the initial point. The virtual mass coef- 
ficient showed little sensitivity to small changes, while damping 
tended to be much more dependent on the point at vrtiich analysis was be- 
gun, Thus, this is a possible cause for the errors being generally 
larger in damping than in virtual mass. 

Because of the use of lAie rotary pitch bearing, friction was neg- 
ligible, In heave a definite friction factor was present. However, as 
discussed in Appendix A , the friction had no effect on the response 
curve and thus virtually no effect on the final values of virtual mass 
and damping. To increase the validity of the results, this friction 
factor should have been included in this analysis program. 



$»2 Character of the Results 



The resultant curves of virtual mass and damping show the general 
form as expected. As predicted by Grim, the curves of virtual mass 
should drop sharply from Infinity to a minimum and then rise to an as- 

symptotic value. The infinite point should occur at approximately a 

(m 

radian frequency of — 9 . = ^ • While the range of frequencies investi- 

g 

gated do not show this point, they do definitely show this trend as it 
is seen that the curves tend toward infinity at frequencies inversely 
proportional to speed, the higher speeds rising at lower frequencies. 
The values of damping also, showed the same trend as for other 
ships. Theory predicts that damping should go to zero at higher fre- 
quencies, However, experimental techniques have shown damping to id.se 
at higher frequencies. This indicates that the theory is only good at 
low frequencies. In some cases the damping appeared to still be 
falling to zero in these experiments. However, the higher speeds 
showed a definite rise. The lower speeds were not tested at higher 
frequencies, as this is the general region of towing carriage equip- 
ment natural frequencies. In this region the validity of the results 
was somewhat questionable an 3 rway. 

For the curves in general a study of the results shows the agree- 
ment of the results in the vicinity of the natural frequency of the 
model to be very good. At frequencies removed from this value the er- 
rors went up considerably. This indicated that this method is accurate 
only near the model's natural frequency. Also, as mentioned earlier, 
the errors in virtual mass were much less than for damping, indicating 



that the procedure may be better for determination of added mass than 
for damping. At the higher frequencies this effect became very pro- 
nounced but was also affected by equipment natural frequencies, 

Kerwin and Narita noted a sudden drop in damping at approximately 
eleven radians per second. This was attributed to the natural frequency 
of the heave rod. In this series of tests a much stiffer heave rod was 
used in an effort to extend the range of frequencies. This sharp dip 
appears to have been eliminated^ but there is enough fluctuation at fre- 
quencies above eleven radians per second to render these results of 
little value. Thus this indicates an even stiffer heave rod would be 
desirable. 

To overcome the problem of frequency range limitations, there is 
also the possibility of altering the model’s natural frequency. This 
was investigated by changing the position of the ballast weights to 
change the radius of gyration. In analyzing the effect of this alter- 
ation, the damping did not appear to change but fell somewhat between 
the two other zero speed conditions. The virtual inertia did change as 
expected and is shown in Figure X, 

The results obtained in pitch did not depend on whether the model 
was initially bow up or stem up. In heave the results were all ob- 
tained with the model lifted up. If a method were used to put pressure 
on the model in order to start with it initially down in the water, it 
would probably produce the same results as when it started lifted up. 
This supports the principle of linearity which was assumed in this en- 
tire procedure. 






5*3 Theoretical Results 



5.3«1 Transfomation Coefficients 

The values of the transfomation coefficients a2n+l 'w®re found 
which mapped a cylinder into the sections of the Mariner, The agree- 
ment between the mapped sections and the actual sections was very good 
in the vicinity of amidships, while the greatest excursions were near 
the ends. For example, at station nine the maximvim difference between 
the mapped section and the actual section was only 0,08 feet with an 
average excursion of 0,039 feet. The worst fit occured at station one, 
which had a maximvim excursion of 0,9315 feet and an average excursion of 
0,1*22 feet# The average excursion for the entire ship was 0,168 feet, 

A jxLot showing the difference in the mapped and actual fom at station 
one is given in Appendix C, 

5,3«2 Comparison of Theoretical and Experimental Results 

From Figures XHI and XIV one can see that the theoretical values of 
virtual mass and damping in heave compare quite favorably with the ex- 
perimentally derived values. The virtual mass «, values show closer agree- 
ment than do the damping values. However, the difference between the 
theoretical and experimental damping is of the same order as in the two 
experimental values found from essentially identical response curves. 

This shows the extreme sensitivity of damping to any minor differences, 
and perhaps strip theory is not a precise enough method to evaluate 
damping, unless very fine station spacing were used. 

In pitch the agreement between the experimental and theoretical 
values of virtual inertia was better than the agreement in damping, 
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However, in both cases the variations were greater than for the corres- 
ponding values in heave. This variation was believed to be the result of 
using strip theory on the Mariner htillform. The Mariner has a bulbous 
bow and a wide flared stem. These sections produced a comparatively 
large value for both added mass and damping. To use strip theory in 
pitch, each section is multiplied by the distance from the center of 
gravity squared and then integrated over the length of the model. These 
extreme bow and stem sections h“ad a very large arm to the center of 
gravity and thus produced very great factors. In the actual case there 
are two reasons why the strip theory is questioned. Near the bow and 
stem the flow is not two-dimensional, as it is near amidship, and thus 
the use of strip theory is not as valid. For example, near the ends the 
water can vent itself not only at the sides, but also at the ends, vrfiich 
should reduce the values found there. These factors may have caused the 
theoretical values to be higher than the values found experimentally. 

As in heave, it is believed that much finer spacing of sections 
should be used, especially near the bow and stem in order to obtain bet- 
ter results for the theoretical values. This will not correct the error 
res\ilting from three dimensional effects, but would reduce errors caused 
by too large an integrating interval in the vicinity of these critical 
regions. 
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Chapter VI 
CONCLUSIONS 



6»1 Eval\iation of Step tfethod 

It is considered from the results of this work that the step dis- 
turbance method is a good method for obtaining added mass and damping 
coefficients. It has several advantages in that the experiment is very 
easy to set up. By use of a device, such as a magnet, it requires only 
low potential devices for power requirements, which are readily avail- 
able, The calculations involved woijld be virtually impossible if done 
by hand, but they are readily adaptable to computer operation. One 
problem that may arise with some models is the problem of wall reflec- 
tions in a narrow tank. This was not a factor in the Mariner tests, 
but by orienting the model to the centerline of the tank at various 
angles, this problem can be overcome. 

It is concluded further that the process is better suited for 
added mass than for damping. This is a factor of the experiment and 
the analysis program. Damping variations were much greater because of 
vibrations or not having the exact starting point in the data. From 
the analysis program, it was seen even with exact generated data, damp- 
ing showed variation from the actual values. However, if care is taken 
to get minimum vibrations and if the data is sampled at fine enough in- 
tervals, it is believed that the process is valid for damping as well 
as for added mass, 

A limiting feature of the experiment is the dependence of the re- 
sults on the model’s natural frequency. This, coupled with the 
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interference of the supporting structure's natural frequency, could 
limit the range of reliable results. By using a very stiff heave rod 
and bracing it, its natural frequency could be made high enough to not 
interfere with the results, 

6.2 Experimental Results 

It is obsearved from the experiments that both added mass and 
damping are speed dependent. However, in the vicinity of the natural 
frequency of the model, this dependency is very slight. The speed de- 
pendence was generally greatest at the lower frequencies, 

C 

Virtual mass was affected by the model's natural frequency, while 
damping showed little dependence. Thus it appears possible to extend 
the useful range of frequencies in damping by testing the model at vari- 
ous natural frequencies. For virtual mass calculations the proper 
natural frequency desired must be used for the model, 

6.3 Theoretical Conclusions 

The method of mapping a cylinder into a ship section is very good 
if sufficient coefficients are used. With the exception of highly ir- 
regular sections, almost perfect agreement between mapped and actual sec- 
tions can be obtained. However, the extent of the calculations make it 
almost impossible to perform them by hand, but they are readily adapt- 
able to computer operation. 

The virtual mass and damping coefficients derived from this method 
were believed to be accurate, as evidenced by the results obtained in 
heave. However, to apply strip theory to the results in pitch, it is 
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recommended that finer spacing be used at the bow and stem, especially 
for irregular sections such as are on the Mariner, By this method a 
more accurate representation of the virtual inertia and damping coef- 
ficients at the ends of the model could be obtained, Tlius, when multi- 
plied by the distance to the center of gravity squared, these sections 
would probably not produce the effect of having the bow and stem domi- 
nate in obtaining the coefficients of virtual inertia and damping for 
the model, ^ 
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Chapter VII 
RECOMMENDATIONS 



7*1 Experimental Technique 

The major problem faced in analyzing the data was obtaining the 
data in digital form* The analog to digital program used was very ef- 
ficient, but there were several mistakes made by the operators which, 
if not made, would have greatly aided the obtaining of the digital 
output. 

The fii*st problem was recording data at too slow a speed and not 
being able to record the digitizing signal. The tape recorder has to 
run at a speed of thirty feet per second to record a one-kilocycle-per- 
second square wave properly. 

Each run took approximately twenty- two seconds to digitize. If 
data were recorded with this spacing between runs, the digitizing con- 
sole should be set up to nxn continuously, and there would be no lost 
time in setting the start of each run when the computer was ready for 
more data. 

In setting up the digitizing console, one has to set up an am- 
plifying and an offset section so as to keep the data between plus and 
minus eight volts. It is best to get as wide a spread as possible in 
order to get the maximum resolution of the data. In order to get this 
spread and not have any data be clipped, it is essential that all runs 
be recorded to the same scale. Thus in setting up the recording equip- 
ment, the zero position and maximum deflection position should be made 
the same on the recording equipment for every run. Otherwise reduced 
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scales must be used to keep all data within the allowable range, with 
some loss in accuracy. 

The final point to be made about digitizing concerns the treat- 
ment of' the magnet parting signal. The method employed by the opera- 
tors in this instance produced several problems that could have been 
avoided. The first problem was that in putting the one kilocycle 
square wave on the tape, the instant of magnet parting was missed in 
some cases as the square wave i/as put on the tape too late. These 
runs had to be digitized a second time to get useful data. This was a 
waste of computer time that was necessary to get the maximum amount of 
data analyzed. Finally, in scanning the record of the pairting signal, 
it is possible that the starting point may have been missed by a couple 
thousandths of a second, irtiis same accuracy could have been obtained 
by simpler methods. The best method appears to be to wire the square 
wave generator between the recorder and the magnet such that when the 
magnet parts, the square wave begins recording. This process should 
have negligible time delay. 

The other big problem was that of random vibrations in the data. 
As stiff a supporting member as possible should be used for supporting 
the ship and supporting the magnet. This would not only reduce the 
vibrations, but also would not limit results because of the equipment's 
natural frequency. Also, the model shovild be accelerated up to speed 
as slowly as possible to minimize the vibrations caused by this accel- 
eration. 

To increase the validity of the results found in heave, a device 
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employing rotary joints and stiff members may have the benefit of less 
friction than the heave rod. The inertia of the mechanism may be dif- 
ficult to determine, but it may be a better mechanism for this type of 
experiment, 

7,2 Future Studies 

The method employed here is considered a good one and could be 
applied to further studies in ship motions. It could possibly be used 
to obtain the values of the coupling coefficients between pitch and 
heave. This could be done with a similar arrangement in which the model 
were allowed both types of freedom, and the two motions recorded. The 

method could also be extended to further studies of damping. This could 

) 

be done by measuring the relation between wave amplitude at the side 
and the amplitude of the model’s oscillation in heave. The same linear- 
ity assumptions could be employed, and the results compared with those 
obtained by forced oscillation tests. 
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APPENDIX A 



DISCUSSION OF EQUHteNT FRICTION 

In choosing the components of the towing apparatus that allowed 
pitch and heave motions of the ship, and in choosing the measuring de- 
vices, it was desired to eliminate or reduce friction to a minimum, 

■a 

Before any connection or measuring plans were made, this was the first 
consideration. 

The problem in pitch was not very great. The only moving con- 
nection, which also contained the measuring device, was the pitch 
beaidng. It is a small bearing rotary device and contributed no no- 
ticeable friction. 

In heave the problem was much greater. The first idea was a 
mechanism employing straight bars connected by small bearing pivots. 

While there would be several such pivots, each would have essentially no 
friction, and the net effect would approach frictionless operation. The 
biggest problem was the system complexity required to insure heave motion 
only. Also, there was a problem of determining the exact system inertia 
in motion. 

Because of the above problems, it was decided to investigate the 
use of a heave rod. Previous heave rods employed had a large amount of 
friction, A new one was constructed using a square rod moving on a 
system of small roller bearings. Friction tests vere conducted on it by 
a weight-dropping method. Various weights were attached to the rod, and 
the acceleration was measured as the rod fell. Comparing this acceleration 
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to one gravity unit^ one could determine the friction in the mechanismo 
Furthermore^ the weights were hung at' fifteen and thiirty degrees to the 
vertical to give the effect of the horizontal force induced by the model 
going through the water at various speeds. The results of these tests 
are tabulated in Table X, 

Table X 

HEAVE ROD FRICTION 

Angle of Weight ^ Average Friction 

0° .385 lbs, 

15° .506 lbs, 

30° .822 lbs. 

To gain insight as to the effect of this friction^ a test was con- 
ducted at zero speed with the heave rod removed. The model was attached 
to the towing carriage such that it was free to oscilldte on the linear- 
syn only. It was set to the same amplitude as was used in the heave 
tests. When the model was released^ and the record of the heave com- 
pared to that with the heave rod^ there was no apparent difference in 
the response. Accordingly^ the effect of heave rod friction was neg- 
lected in the final calculations. 

As this test was conducted at zero speed and as the weight drop 
tests showed considerable increase of friction with a horizontal coin- 
ponentj this may have been a source of slight error at higher Froude 
numbers. However^ this effect was considered secondary and was neg- 
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APPENDIX B 



EFFECTS OF WALL REFLECTIONS 

A question raised early in the testing was how valid the results 
would be or whether the effects of wave reflections in a narrow tank 
would be such as to render the response meaningless. Because of this 
possibility, two procedures were available to get better results at 
zero speed. One was to test the model at different orientations in the 
tank, and if this proved unsuccessful, to make the tests in a larger 
medivun, such as a swimming pool. 

The model was first tested parallel to the tank. The heave re- 
sponse obtained is shown in Figure XVII, The model was then tested at 
an orientation perpendicular to the tank, and the response is shown in 
Figure XVIII, Similar results were obtained for pitch. In both cases it 
is noted that the model oscillations damped to zero before any wall re- 
flections affected the motion of the model. It is true that the ex- 
tremely long waves returned to the model almost instantly, but their 
effect is negligible, and it is the wave of the same order of magnitude 
as the ship length with which we are primarily concerned. 

It is interesting to note that there is a considerably longer 
zero region when the model was placed across the tank than when it was 
oriented parallel to the tank. This indicates that wall reflections 
must be considered in this type of analysis, but fortunately, for the 
Mariner model the oscillations completely damped out prior to wall re- 
flections, Thus the results of the analysis were essentially the same 
for each condition. 
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APPENDIX C 



DETERMINATION OF COEFFICIENTS OF TRANSFORMATION 

The determination of the coefficients of transformation of the 
semicircle in one complex plane to the ship section in the second com- 
plex plane pixived to be very adaptable to computer operation. The IBM 
7091; computer at the Massachusetts Institute of Technology was used 
for this purpose in this project. In general, the mapped cylinder gave 
a very close approximation to the actual ship section. The method used 
consisted of selecting tentative values of a]^, a^, and a^ by reference 
to Landweber and Macagno (9)« The computer program was then entered 
with these values and offsets of the actual ship section. An iteration 
procedure written into the program generated new values of the coef- 
ficients and offsets which were compared with the actual ship section. 
The iterations were continued until a reasonably small least mean 
square error between calculated and actual offsets was obtained. The 
maximum deviation used in this project occured at station one. Fig- 
ure XIX in this appendix compares the actual station one idLth station 
one derived from the transformation coefficients. 
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COMPARISON OF THE ACTUAL MARINER STATION ONE 
WITH THE CALCULATED SECTION 




APPENDIX D 



NEW METHOD FOR DETERMINATION OF TRANSFORMATION COEFFICIENTS 



After the main body of work on this project was completed, a new 
method for determining the transformation coefficients, a2n+l> 
vised. This method, together with the nev; values of the transformation 
coefficients for the stations of the Mariner hullform, are presented in 
this appendix. 

Equation (5) represents the conformal transformation of the unit 
circle to the section in question. In equation (5) 
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X + iy = ire 
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(18) 



In this work the y axis points down, and the angles 9 and 3 are measured 
from the vertical. X is the offset of the actual form, and y is the 
waterline. Therefore, when 9=3= 90° ^ x is the half -beam of the sec- 
tion, and y = 0, When 9 = 3 = 0°, x = 0, and y is the draft of the sec- 
tion. 

By substituting equations (I6) and (17) into equation (5) and sepa- 
rating the real and imaginary parts, the following equations are obtained: 

N 

X = psing + ^ (-1)’^ “§71 + 1)9 (19) 
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y = pcos0 + 



N 



cos(2n + 1)0 



( 20 ) 




n=0 



(- 1 ) 



n+1 



^2n+l 

2n+l 

P 



It is then necessary to solve for 0 and satisfy these equa- 

tions. 

If M points on a section outline are available, there are M pairs 
of equations of the above form, or a total of 2M equations. The values 
of a2n+l sane in all of these equations, but the values of 0 

will be unique for each given point. 

If the half -bean and draft of the section are known, an additional 
equation is available. As mentioned above, at 0 = 90 °, 



half -beam = b 



At e = 0°, 



P + 



V"' 2nH-l 
n=0 ^ 



draft = H 



N 

P + 

n=0 



®2n+l 

^2n-i-l 

P 



( 21 ) 



( 22 ) 



These results lead to the equation 

^ 0 . 

p(l - |) + (1 + |) + (1 - |) + .... =0 ( 23 ) 

P 

■Wte now have 2M+1 equations which must be satisfied. 

In the new method of solving for the coefficients a2n-‘-l^ 
ordinates of the M points, as well as the half -beam and draft of the 



section, are inputs to the computer program* Three alternatives are now 
available. If estimated values of the coefficients 32^+1 are available, 
they may be entered, In lieu of this, if estimates of the angles 6 are 
available, they may be entered as inputs to the program. If neither the 
coefficients nor the angles 9 are available, the program will generate 
an estimate of the an^es 0. 

Assuming that either the angles 9 are inputs or that values of © 
are internally generated, the program solves the 2M«-1 equations for P 
and a2n*-l ^ give a least mean square error solution between the gener- 
ated section and the actual ship section. The values of the coeffic- 
ients are then used with input values of x and y to determine a new 
value of 9 for each point. The program then re-enters the least mean 
square process and iterates until the desired number of iterations is 
completed. If values of ^2n*l inputs, the program computes values 
of 9 before entering the mean square process. 

It should be noted that the program requires the number of co- 
efficients desired and the number of iterations desired as inputs. 

As stated above, the program finds values of p. However, since 
we are interested in the case of the unit circle, the program is nor- 
malized with p equal to one. 

Using this new method, new values of were determined for the 

Mariner huUform, These coefficients are presented in Tables XI and XII. 

It is believed that this improved approximation to the hullform 
will not appreciably affect the values of damping and virtual mass for 
the Mariner, However, in future work it is recommended that the new 
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method presented above be used, 

A listing of the analysis program in MAD computer language is at- 
tached to this appendix. 

Table H 

TRANSFORMATION COEFFICIENTS 



Station 






0 


-.858571 


-.118183 


1 


-.675965 


-.085363 


2 


-.U79528 


-.Oii3086 


3 


-.286588 


-,022750- 


U 


-.117882 


-.0208iiii 




-.OOOU35 


-.035251 


6 


+.076025 


-.053853 


7 


+.121668 


-.078578 


8 


+ .li;399h 


-.102216 


9 


+.150007 


-.1190i;0 


10 


+.150031 


-.12i;7ii6 


11 


+.l5Hi;li 


-.122959 


12 


+.150825 


-.110575 


13 


+.lh8566 


-.085017 


11; 


+,139671 


-.OU8057 


15 


+.117lil0 


-.OO2I430 


16 


+,069663 


+.053825 


17 


-.0269U9 


+.106876 


18 


-.216518 


+ ,lli6021 


19 


-.562587 


+,131206 


20 


+.227269 


+.251236 



!i 






-,018627 


+.006032 


-.000315 


-.51353I4 


+.006935 


+.005787 


+.OIOI33 


+.003225 


+.013838 


+ .OI9I+79 


+.003222 


+.016652 


+.019821 


+.005278 


+.010582 


+.017899 


+.002628 


+ .00781iU 


+.OI6527 


+ .001i;57 


+.005735 


+ .0107l;7 


-.OOO67I; 


+.002915 


+ .00ia09 


-.001911 




-.001388 


-.002853 




-.OOi;275 


-.00i;280 


+.001686 


-.003135 


-.ooiaoii 




-.001117 


-.001822 




+.005537 


-.OOi;659 




+.016075 


-.OOU566 


+.001582 


+ .02ii3ii5 


-.001183 


+.006575 


+.032166 


+ .0000li2 


+.00933U 


+.0ii35l2 


-.001900 


+ , 006618 


+ .07iai5 


+,012576 


-.003023 


+ .0982li5 


+.0^7890 


+.016198 


-.12U012 


+ .Oli8271 


-.109952 
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0 



•I 



1 

2 

3 

h 

5 

6 

7 

8 

9 

XO 

11 

12 

13 

lU 

15 

16 

17 

18 

19 

20 



Table XU 



COMPARISON OF ERRORS WITH THOSE FOUND BY PROCESS 
U§ED TiJ THE MAIN BODY OF TfiE PROJECT 



New Method 



Main Body- 



Average 

Excursion 



Maximum 

Excursion 



Average 

Excursion 



Maximtun 

Excursion 



0.07 


0.11 






0.09 


0.17 


0.32 


0.93 


0.16 


0.33 


0.17 


0.32 


0.13 


0.25 


0.20 


0.51; 


0.12 


' 0.23 


0.11 


0. i ;0 


0.10 


0.22 


0.13 


0.32 


0.08 


0.17 


0.08 


0.17 


0.06 


0.13 


0.06 


0.18 


0.03 


0.07 


0.07 


0.13 


0.02 


0.07 


0.02 


0.08 


0.05 


0.12 


0.08 


0.23 


0.06 


0.11 


0.06 


0.11 


0.03 


0.09 


0.03 


0.09 


0.03 


0.05 


o.ou 


0.08 


0.03 


0.09 


o.oU 


0.10 


0.09 


0.25 


0.19 


0.1;0 


0.08 


0.19 


0.2 i | 


0.55 


0.12 


0.27 


0.19 


0.52 


0.06 


0.11 


0.22 


0.1;6 


0.15 


0.30 


0.29 


0,69 


0.07 


0.13 


0.07 


0.20 
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RCOEFFICIENT DETERMINATION FOR HULL FORMS 
DIMENSION B(62S.BDIM) ,THETA(62S,TDIM) »X(30) »Y(30) »ABAR(23 ) ♦ 

1 APR I (15)»BETA(30).CRlT(20).A(15)»WYE(30).tiX(30)»tRRX(30)t 

2 ERRY(30) »NAME( 12) » DEGR E E ( 62 t? * OD I M ) .ALPHA (30) » DUM ( 2 2 .U I MD ) 
INTEGER I .J.K.M.N.Qt COUNT .PWR. SIGNAL 

VECTOR VALUES BDIM = 2.1.25 

VECTOR VALUES TDIM = 2.1.25 

VECTOR VALUES DDIM = 2.1.25 


START 


VECTOR VALUES DIMD = 2.1.15 

READ FORMAT T I T L 1 . NAME ( 1 ) . . . NAME ( 1 2 ) 

VECTOR VALUES TITLl = $(12C6)*$ 

PRINT FORMAT T I TL2 . NAME ( 1 ) . . . NAME ( 1 2 ) 

VECTOR VALUES TITL2 = $ ( IHl . 12C6 ) *S 

READ format INPUT. H ♦ m , N .0 . BE AM . S I GN AL 

VECTOR VALUES INPUT = i ( F 1 0 . 6 , 3 I 2 » F 1 0 . 6 . I 2 ) 

BDIM(2) = N+2 

DIMD(2) = N+1 

TDIM(2) = 0+1 

D0IM(2) = 0+1 

READ FORMAT OFFSET, X(1),,.X(M) 

VECTOR VALUES OFFSET = f(5F10.4)*$ 

READ FORMAT WL . Y(1)..,Y(M) 

VECTOR VALUES WL = S(5F10.4)*S 

KAY = BEAM/H 

El = l.OE-06 

WHENEVER SIGNAL #E, 0 

READ FORMAT ANGLE. EX ( 1 ). . . E X ( M ) 

THROUGH BEEl. FOR' I = 1.1. I.G.M 


BEEl 


DEGREE! I .1 ) = EX( I ) 

VECTOR VALUES ANGLE = S(5F10.5)*$ 
THROUGH TWENTY. FOR I = l.l.I.G.M 


TWENTY 


THETA(I.l) = DEGREE! I .1 ) /57. 295780 
TRANSFER TO TWELVE 
OR WHENEVER SIGNAL .E. 1 
READ FORMAT AV/L. A!1)*..A(N) 

VECTOR VALUES AVAL = $ ! S4 ♦ 4E 1 4 ,6 ) *$ 


thirty 


THROUGH thirty. FOR I = 1 . 2 . I .G. ! 2*N-1 ) 
A! 2*1-1 ) = A! I ) 

DEN = 1. + A!1 ) 

THROUGH FORTYf FOR I = 3 .2 ♦ I .G. ! 2*N-1 ) 


FORTY 


DEN = DEN + A ! I ) 
RO = BEAM/DEN 


FIFTY 


THROUGH FIFTY, FOR I- = 1 .2 . I .G, ! 2*N-1 ) 
APRI ! I ) = A! I )*R0 
TRANSFER TO TWELVE 
OTHERWISE 
CONTINUE 

END OE CONDITIONAL 
R 
R 

THROUGH INITL. FOR I = l.l.I.G.M 


INITL 

TWELVE 


THETA!I,1) = ATAN. (X( I ) /Y! I ) ) 

WHENEVER Y!I) .LE. El, THETA!I.1) = 1.57080 
THROUGH ONE. EOR J = l.l.J.G.Q 
WHENEVER SIGNAL .E. 1 .AND. J .E. 1. TRANSEER TO B1 



1, TRANfFER TO B1 



1 






TWO 

Three 



FOUR 



FIVE 



SIXPRI 
S I X 



EIGHT 






Q 

P 

THROUGH TWOt FOR I = !,!♦ I.G.M 
B ( I * 1 ) = 51 N. ( theta ( I , J ) ) 

B( I .2) = B( I .1 ) 

B ( I ♦ ( N + 2 ) ) = X ( I ) 

B((I+M).l) = COS. (THETA( I tJ) ) 

B( ( I+M) .2 ) = -F ( ( I+M) , 1 ) 

B( ( I+M) . (N+2) ) = Y( I ) 

THROUGH TWO. FOR X = 3.1. K.G.(N+1) 

B(I.K) = ( (-1 . ) .P.K )*SIN. ( ( ( 2*K)-3 )*THETA( I . J ) ) 

B((I+M).K) = - ( ( -1 , ) .P.K ) *COS. ( ( ( 2*K ) -3 ) *ThETA( I . J ) ) 
through three, for K = 1 .1 .K.G. ( N+1 ) 

B((2*M+1).K) = 1, + ( (-1 . ) .P.K)*KAY 

B ( ( 2*M + 1 ) . (N + 2 ) ) = 0, 

R 

R 

EXECUTE GLSOH. ( B.ABAR .DUM. ( N+1 ) . ( 2*M+1 ) ) 

R 

PRINT RESULTS ABAR ( 1 ) . . . ABAR ( N+1 ) 

R 

RO = ABAR(l) 

through four, for I = 2.1. I.G. (N+1) 

APRI ( (2*1 )-3) = AB^R( I ) 

R 

R 

THROUGH FIVE. FOR I = 1.2. I . G. ( ( 2 *N ) - 1 ) 

A ( I ) = APR 1(1) /RO 
R 
R 

THROUGH SIX. FOR I = 1.1. I.G.M 
EX(I) = (RO + APRI ( 1 ) )*SIN. ( THETA( I . J ) ) 

PWR = 0 

DEGREE(I.J) = 57.2958*THETA( I .J) 

WYE(I) = (RO-APRI ( 1 ) )*COS. ( THETA( I .J) ) 

THROUGH SIXPRI. FOR K = 3.2. K. G. ( ( 2*N ) -1 ) 

PWR = PWR + 1 

WYE(I) = WYE(I) - ( (-1 , ) .P.PWR )*APR I ( K )*COS. ( K + THETA ( I . J ) ) 
EX(I) = EX(n + ( (-1. ) .P.PWR )*APRI (K)*SIN. (K*THETA ( 1 .J) ) 
FRRY(I) = .APS.(WYE(I) - Y(I)) 

ERRX ( I ) = .ABS. ( EX ( I )-X ( I ) ) 

R 

R 

ALPHA(O) = 0. 

BETA(O) = 0. 

THROUGH EIGHT. FOR I = 1.1. I.G.M 
ALPHA(I) = ALPI-'A(I-l) + ERR X ( I ) *ERRX ( I ) 

BETA(I) = BETA(I-l) + ERRY ( I )*ERRY ( I ) 

R 

R 

CRIT(J) = SORT, ( ALPHA(M) + BETA(M)) 

EXECUTE OUTPT, ( J . RO . A . CR I T . X . EX . ERRX . Y . WYE . ERR Y . DEGREE . M . N ) 
whenever J.E.O. transfer to FINIS 
R 
R 



I 



>k. 





through CEEl* FOR I = 1»2» I.G.(2*N-1) 

PRINT FORMAT ABC. A( T ) 

VECTOR VALUES ABC = S(1H ,E14.6)*S 
PRINT FORMAT LMS. CRIT(J) 

VECTOR VALUES LMS = S(21H LEAST MEAN SQUARE = .F14,6)*S 

PRINT COMMENT $0 MEASURED X(I) CALCULATED X(I) 

1 ERROR measured Y(I) CALCULATED Y(I) Y(I) 

20R THETA(I)$ 

THROUGH TEN. FOR I = 1.1. I.G.M 
PRINT FORMAT OFVAL ♦ X ( I ) ♦ EX ( I ) , ERRX ( I ) . Y ( I ) . WYE ( I ) . ERRY ( 

1 DEGREE(I.J) 

VECTOR VALUES OFVAL = $(1H ♦ F 14 . 6 . S5 ♦ F 1 4 . 6 ♦ S 5 . F 14 . 6 . S5 . 

1 F14.6.S5,F14,6.S5.F14.6,S5.F14.6)*S 

FUNCTION RETURN 
END OF FUNCTION 

external function (A.X.B.N.M) 

entry to GLSOH. 

INTEGER I.J.K.N.M 

THROUGH QQl. FOR I = 1.1. I#G*N 

X ( I ) = 0. 

THROUGH QQl. FOR J = 1.1» J.G.M 
X(I) = X(I) + A( J.I )*A( J. (N + 1 ) ) 
through QQ2. FOR I = 1.1. I.G.N 
THROUGH Q02. FOR J ='1.1. J.G.N 
B( I .J) = 0. 

THROUGH Q02. FOR K = 1.1. K.G.M 
B(I.J) = B(I»J) .+ A(K.I)*A(K.J) 

THROUGH Q03. FOR I = 1.1. I.G.N 
WHENEVER B( I . I ) .E. 0. 

THROUGH QQ4. FOR K = (I + D.l. K.G.N 
THROUGH QQ5. FCR J = I.l. J.G.N 
TEST = B( I »J) 

B ( I ♦ J ) = B (K ♦ J ) 

B(K.J) = TFST 
TEST = X ( I ) 

X( I ) = X(K) 

X(K) = TEST 

WHENEVER B(I.I) .NE. 0.. TRANSFER TO QQ6 
END OF CONDITIONAL 
STORE = B ( I ♦ I ) 

THROUGH QQ7. FOR J = I.l. J.G.N 

B(I.J) = B(I.J) /STORE 

X( I ) = X( I )/STORE 

THROUGH QQ3, FOR K = 1.1. K.G.N 

WHENEVER K . E . I . OR . B ( K . I ) . E . 0 . ♦ TRANSFER TO QQ3 
STORE = B (K. I ) 

THROUGH QQ8. FOR J = I.l,. J.G.N 
B(K,J) = B(I.J)*STORE - B(K.J) 
whenever K.L.I. B(K.J) = -B(K.J) 

X(K) = X( I )*STORE - X(K) 

WHENEVER K.L.I. X(K) = -X(K) 

CONTINUE 
FUNCTION RETURN 
END OF FUNCTION 



X ( I ) 
ERR 

I ) . 



^iinin 3 Mariner class hullfor 
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